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ABSTRACT 

We  find  that  changes  in  threshold  voltage  induced  by  negative  bias  temperature  stressing  of  p- 
channel  field  effect  transistors  with  HfSiON  gate  dielectrics  are  modulated  by  the  drain  voltage, 
in  measurements  wherein  the  drain  current  is  measured  during  stressing.  This  effect  is  not 
observed  in  SiC>2  gate  devices.  Short  channel  effects  are  excluded  as  explanations,  leading  us  to 
conclude  that  positive  charge  in  the  dielectric  stack  is  laterally  mobile  and  is  conducted  out  of  the 
insulator  via  the  drain.  Further,  a  simple  qualitative  model  of  charging  kinetics  allows  us  to 
extract  the  density  of  interface  states  as  a  function  of  time,  and  shows  that  these  defect  densities 
of  the  order  of  10  cm'  after  hundreds  of  seconds.  These  values  are  consistent  with  observations 
on  pure  SiC>2  gate  insulator  devices. 
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1.  Introduction 


Accurate  determination  of  the  susceptibility  of  metal  oxide  semiconductor  field  effect 
transistors  (MOSFETs)  to  degradation  via  the  Negative  Bias-Temperature  Instability  (NBTI)  is 
of  crucial  importance  for  future  generations  of  Si  based  devices,  since  this  has  been  recognized 
as  one  of  the  critical  reliability  issues  in  advanced  CMOS  devices  [1].  Complicating  this  task  is 
the  fact  that  there  appears  to  be  no  established  generic  measurement  protocol;  indeed,  they  vary 
considerably  from  one  research  group  to  another,  making  cross-referencing  of  research  results 
difficult.  In  our  own  recent  investigations  [2],  to  minimize  errors  resulting  from  trapped  charge 
relaxation  during  the  MOSFET  characteristic  acquisition  sequence,  we  have  adopted  the  practice 
of  measuring  changes  in  the  drain  current  Ids  whilst  maintaining  a  constant  bias  stress  level  Vgs 
on  the  gate  electrode.  The  device  threshold  voltage  shift  AVti,  under  this  NBTI  stress  can  then  be 

Af ds 

determined  from  the  fractional  current  variation  Alds/[ds  fds  given  certain  approximations 
which  we  have  discussed  at  some  length  [2]. 

Under  usual  conditions,  Ids  can  be  determined  in  either  the  linear  or  saturation  modes  of 
operation  [3]  of  the  MOSFET.  Following  reference  [3],  Ids  can  be  approximated  by 


lds  =(W/L)vaffCox[(vgs-Vth)vds-\vl\,  ld. 

where  W  and  L  are  the  device  inversion  channel  width  and  length,  pcn  is  the  effective  inversion 
channel  carrier  mobility,  and  Cox  the  gate  insulator  capacitance.  The  drain-source  voltage  VdS 
biases  the  device  in  the  linear  regime  when  '7<fr  »  0 rgs  —  )  |FaJ  -  Vth  |  and  in  the  saturation 


regime  when  |*  d*  -  (K?s  l  -  PtJ.  Assuming  that  some  physical  phenomenon  (such  as 

NBTI)  induces  a  given  AVti,  and  an  associated  Aids,  then  from  Eq.  (1)  it  is  easy  to  show  that 
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so  that  for  a  maximum  measurable  effect  it  is  more  sensible  to  measure  in  the  saturation  regime. 
However  in  practice,  measurements  are  usually  made  in  the  linear  regime  where  VdS  is  small. 
Logic  for  the  latter  may  be  the  possibility  of  non-uniform  charge  generation  or  trapping  in  the 
gate  dielectric  or  at  the  Si/dielectric  interface  resulting  from  a  non-uniform  vertical  (z  direction) 
electric  field  Ez(y),  when  VdS  is  not  small  and  so  the  semiconductor  surface  potential  becomes  a 
function  of  distance  along  the  inversion  channel  (y  direction). 


In  the  work  reported  here  we  have  performed  NBTI  measurements  on  Si02  and  HfSiON 
gate  dielectric  p-channel  MOSFETs  while  stressing  in  either  the  linear  or  saturation  drain  current 
regime.  We  have  observed  what  appears  to  be  an  important  measurement  protocol  issue,  at  least 
in  the  HfSiON  devices  we  have  studied:  an  unexpected  dependence  of  the  NBTl-induced 
threshold  voltage  shift  on  the  drain  voltage  applied  during  stress.  We  argue  that  this  is  not  due  to 
lateral  non-uniformities.  That  different  VdS  could  affect  NBTI  stress  measurements  is  important 
to  note  both  from  the  point  of  view  of  physical  understanding  of  the  effect  and  from  the  point  of 
view  of  deriving  predictive  formulae  for  reliability  lifetimes.  Furthermore,  based  upon  our 
observations  of  the  effect  of  VdS,  we  propose  a  simple  qualitative  model  of  charging  kinetics  that, 
while  approximate,  allows  us  to  also  extract  the  NBTI-induced  interface  state  density  as  a 
function  of  time. 


2.  Experimental  Procedure 

The  high-k  samples  used  in  this  investigation  were  pMOSFETs  in  which  the  dielectric 
stack  consisted  of  a  nominal  1  mn  of  SiC>2  onto  which  was  deposited  2  nm  of  HfSiO  using 
atomic  layer  deposition.  The  stack  was  subsequently  nitrided  and  the  resultant  atomic 
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composition  estimated  to  be  SiCb  (70%),  HfC>2  (30%)  and  N  (15%)  with  respect  to  the 
concentration  of  starting  HfSiO  material.  The  effective  dielectric  constant  for  the  stack  was 
approximately  13.  A  detailed  description  of  the  gate  dielectric  manufacturing  process  is  provided 
in  reference  [4].  Devices  with  gate  lengths  of  1.0  pm  and  0.7  pm,  both  with  0.5  pm  gate  widths, 
were  measured;  the  two  gate  length  sets  gave  similar  results.  At  least  20  HfSiON  FETs  were 
used  in  our  studies  to  ensure  reproducibility  of  the  observed  phenomena.  Additionally,  p-channel 
devices  with  2.6  nm  thick  SiC>2  gate  dielectric  were  fabricated  on  silicon-on-insulator  substrates 
as  described  elsewhere  [5].  The  measured  devices  had  typical  channel  lengths  ~  0.5  pm  and 
widths  <  2  pm. 

Drain  current  measurements  were  made  using  a  Keithley  4200  SCS  measurement  system 
in  “continuous”  mode  [2]  at  a  gate  voltage  of  -2  V,  and  an  arbitrarily  chosen  temperature  of  175 
°C  .  Because  the  NBTI  is  generally  observed  to  be  smaller  in  pure  SiC>2  gate  dielectric  devices, 
these  were  measured  at  195  °C  to  enhance  the  effect.  The  drain  voltage  V<js  was  varied  over  the 
range  -0.1  V  to  -2.0  V,  but  was  constant  during  measurement  on  a  particular  device.  Source, 
substrate  and  body  were  grounded.  In  the  continuous  mode,  the  drain  current  IdS(t)  was  measured 
periodically  (approximately  every  0.75  s)  without  interrupting  the  stress.  This  nominally  avoids 
the  likelihood  of  trapped  charge  relaxation  [6]  which  is  commonly  observed  in  experiments 
where  the  stressing  voltage  is  removed  in  order  to  perform  an  IdS  versus  Vgs  measurement.  From 
the  drain  current  data,  the  threshold  voltage  shift  AVti,(t)  was  extracted  using  the  simple  square 
law  model  (Eq.  1) ,  which  yields 


(3) 


where  'V*  is  the  threshold  voltage  prior  to  NBTI  stress  (~  -0.5  V  at  175  °C  for  the  HfSiON 


devices  and  +1.3  V  at  195  °C  for  the  SiC>2  devices),  is  the  drain  current  at  t  =  0,  and  Vgs  and 
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Vds  are  the  gate  and  drain  voltages.  Eq.  (3)  holds  provided  ^ds  »  0  W  “  ) 

|I^,.  -  Vth  | 'as  >  —  ^'tho ;  for  more  negative  VdS,  _  ^  :h )  l$*  —  ‘the  must  be  substituted  for 

VdS.  Note  that  for  the  measurements  on  the  SiC>2  with  large,  positive  Vti„  Vds  =  Vgs  =  -2  V  does 
not  strictly  meet  the  saturation  definition  condition.  However,  we  performed  Ids(Vds) 
measurements  for  Vgs  =  -  2  V  and  observed  that  the  current  was  within  4%  of  the  saturation 
value. 

The  protocol  described  here  is  designed  to  minimize  charge  loss  during  measurement. 
The  largest  uncertainty  in  these  measurements  is  in  determining  the  initial  source-drain  current 
(Ids°)  since  the  system  takes  1  -  2  seconds  to  generate  the  first  data  point.  We  estimate  the 
uncertainty  in  lds°  to  be  of  the  order  of  0.1%,  causing  errors  of  a  few  percent  in  the  extracted 
AVti,.  We  also  ignore  any  changes  in  mobility,  which  we  have  previously  argued  to  be  a 
reasonable  approximation  for  these  samples  [2]. 


3.  Results  and  Discussion 

Normalized  change  in  drain  current  l  ds  Alds(t)/lds  measured  during  stressing,  and 
the  resultant  threshold  voltage  shifts,  are  shown  for  HfSiON  gate  devices  in  Fig.  1 .  All  samples 
were  stressed  at  175  °C  and  Vgs  =  -  2.0  V,  and  at  varying  Vds.  As  indicated  in  Eq.  (2),  the  square 
law  model  predicts  that  for  a  given  AVu„  the  normalized  change  in  drain  current  Aids/  Ids° 
measured  in  the  saturation  region  <  !'j5j  -  1  'itho)  should  be  twice  as  big  as  in  the  linear 
region  (VdS  ~  0).  In  any  event,  AV,ii  should  not  depend  on  Vds.  What  we  observe,  however,  is  that 
in  saturation  (Vds  <  -1.5  V)  the  change  in  drain  current  is  in  fact  much  less  than  in  the  linear 


region  (Vds  =  -0. 1  V).  As  a  result,  AVm  is  found  to  be  smaller  by  a  factor  of  more  than  4  when 
NBTI  stress  is  performed  in  saturation  as  compared  with  the  value  estimated  from  the  linear 
region  data  (Fig.  lb).  Intermediate  values  of  VdS  have  a  corresponding  effect  on  AVti,. 

We  have  also  performed  NBTI  measurements  at  1 95°C  on  p-channel  MOSFETs 
constructed  on  silicon-on-insulator  substrates.  Ias(t)  measurements  were  made  with  the  source- 
drain  voltage  biased  in  the  linear  regime  ( -0.2  V)  and  in  the  saturation  regime  (-  2V).  The 
results  for  AWIds°  are  shown  in  Fig.  2a.  and  after  conversion  to  AVt|,(t),  in  Fig.  2b.  As  argued 
above,  we  observe  that  AIdS/  IdS°  in  the  saturation  regime  is  ~  2  x  the  variation  in  the  linear 
regime.  This  result  is  in  agreement  with  what  is  anticipated  in  a  p-channel  MOSFET  but  is  in 
complete  contrast  with  what  we  have  observed  for  HfSiON  gate  devices.  This  conclusion  for 
Si02  based  devices  is  further  supported  in  other  published  data  [6,7]. 

3.1.  Short  channel  effects 

The  first  explanation  of  the  drain  voltage  dependence  of  AVu,  which  comes  to  mind  is  a 
short-channel  effect  in  which  charge  is  created  non-uniformly  at  higher  drain  voltages.  If  that 
were  the  case,  one  would  expect  to  observe  it  in  SiC>2  devices  as  well,  but  we  do  not.  There  are 
other  reasons  why  non-uniform  charging  is  not  likely  to  explain  the  observations.  We  have 
performed  measurements  of  IdS  as  a  function  of  Vgs  in  the  sub-threshold  regime  for  our  devices. 
Sub-threshold  current  plots  of  our  1  pm  devices  measured  with  VdS  varying  from  -0. 1  to  -2.0  V 
show  no  measureable  shift  in  Vti,°.  This  is  in  fact  to  be  expected  since  the  doping  density  in  the 
devices  was  ~10  cm'  and  the  equivalent  oxide  thickness  was  approximately  2  nm.  Thus  the 
devices  are  very  “long  channel-like”,  so  drain  voltage  effects  are  confined  to  a  region  very  close 
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to  the  drain.  Furthermore,  we  note  that  there  is  a  sizeable  change  (Fig.  lb)  in  AV,h  even  for  lower 
drain  voltage  (-0. 1  V  >  V<]S  >-1.0  V),  i.e.,  well  before  saturation. 

We  can  place  the  “long  channel-like”  argument  on  more  sound  physical  footing  using  a 
quasi-two-dimensional  model  [8].  The  model  allows  calculation  of  the  surface  potential  along 
the  channel  in  the  Si  substrate  Vs  while  accounting  for  2D  effects  present  near  the  source/drain 
contacts  in  small  devices.  Specifically,  we  have  [8] 

dmaxdExb 0  .  _  Yas~  vfb  -  V,(y) 


sSi' 


+  S„ 


=  qNBd, 


(4) 


V  dy  tox 

The  z-  and  y-directions  here  correspond  to  the  vertical  and  lateral  directions,  so  that  Ez(y)  is  the 
vertical  oxide  field,  which  is  a  function  of  distance  along  the  channel.  The  dielectric  permittivity 
in  the  silicon  and  in  the  oxide  are  esi  and  eox,  Vgs  and  Vfb  are  the  applied  gate  voltage  and 
flatband  voltage  [3],  Nb  and  tox  are  the  substrate  doping  density  and  oxide  thickness,  d,nax  is  the 
maximum  depletion  depth  in  the  substrate  due  to  the  drain  voltage  (and  is  given  by  the  depletion 
approximation),  and  r)  is  a  fitting  parameter,  which  we  set  to  1 .  The  solution  for  the  surface 
potential,  with  boundary  conditions  Vs(0)  =  Vbi  and  VS(L)  =  Vas  +  Vbi  is 

fL  -) 


V. 


sinhi j)  sinhC — j-M 

SO’)  =  vA  +  (yhi  +  vas  -  vsL) — Hj-  +- <y¥i  -  v}L)  — 

sinhy jj  siJih  (yj 


(5) 


Here,  —  ~  !  tho  +  V*  ^ sL  -  ~  z< Vb  where  Kta  is  the  long-channel  threshold 

voltage,  <pB  is  the  Fermi  potential  in  the  silicon  bulk,  and  Vgs  and  Vbi  are  the  applied  gate  voltage 


and  built-in  potential.  The  length  parameter  ^  is  given  by  /  = 


l  = 


fajtoxd 


max 


<  soxV 


Zox>l 


We  have  applied  the  model  to  the  geometry  of  our  devices,  using  Vbi  =  I  V,  Vfb  =  0,  and  I'Vi  =  - 
0.5  V,  the  oxide  permittivity  appropriate  to  SiC>2  with  an  oxide  thickness  tox  of  2  nm,  and 
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semiconductor  parameters  at  1 75  °C.  From  the  surface  potential  we  have  calculated  the  oxide 

Vas  -  l'5(v) 

E2(y)  =  — - — 

field  as  a  function  of  distance  along  the  chamiel  as  tox  .  A  plot  of  Ez(y)  is  shown 

in  Fig.  3a  for  the  case  of  VdS  =  -0.1  V  and  -2  V.  We  observe  that  Ez(y)  is  very  uniform  along  the 
channel  for  both  cases  and  that  for  the  -2V  case  significant  variation  occurs  only  in  the  last  0.05 
pm.  Therefore,  whether  in  the  linear  regime  or  the  saturation  regime,  our  devices  should  have 
experienced  a  uniform  Ez(y)  profile  over  nearly  the  entire  channel. 


3.2.  Non-uniform  trapped  charge  generation 

The  previous  estimates  reveal  that  the  vertical  electric  field  is  essentially  uniform  along 
the  channel  length  in  both  linear  and  saturation  cases.  However,  we  cannot  a  priori  reject  the  idea 
that  the  trapped  charge  generation  mechanism  itself  in  some  way  “amplifies”  the  non-uniformity 
of  the  electric  field  and  thereby  gives  rise  to  substantially  less  charge  in  one  case  as  opposed  to 
the  other.  Various  authors  have  endeavored  to  estimate  the  vertical  electric  field  dependence  of 
the  NBT1  charge  generation  rate  G(EZ)  resulting  from  the  NBTI  phenomenon.  We  arbitrarily 

chose  one  which  has  the  form  [9]  & (Ez)~Eze]  where  y  has  a  value  of  0.6.  We  show  in  Fig.  3b 
a  plot  of  the  generation  rate  as  a  function  of  distance  along  the  device  channel  assuming  the  Ez(y) 
variation  also  shown  in  Fig.  3a  -  the  calculation  assumed  Vds  =  -  2V.  The  generation  rate  of 
trapped  charge  due  to  NBTI  is  again  found  to  be  relatively  uniform  along  the  channel.  We  are 
therefore  led  to  conclude  that  short  channel  effects  in  general,  and  non-uniform  charging  in 
particular,  cannot  give  rise  to  a  level  of  trapped  charge  generation  substantially  different  from 
that  one  would  anticipate  if  Ez(y)  were  entirely  uniform  as  in  the  case  of  linear  regime  device 
operation  (low  Vds). 
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3.3.  An  alternative  explanation 

To  recap  the  situation,  we  have  performed  NBTI  measurements  on  devices  operated  in  the  linear 
and  saturation  regimes.  Surprisingly,  the  experimental  data  on  high-k  based  devices  gives 
evidence  that  charge  trapping  is  substantially  larger  in  the  former  operational  mode  than  in  the 
latter.  Having  eliminated  short  channel  and  non-uniform  charging  effects,  we  look  for  an 
explanation  in  the  nature  of  the  dielectric  itself.  Noting  that  high-K  dielectrics  tend  to  trap  charge 
more  readily  than  Si02,  we  postulate  the  following:  NBTI  induces  a  high  density  of  positive 
charge  in  the  dielectric  near  the  interface  with  the  Si  substrate,  perhaps  in  the  thin  “SiCV’  layer 
itself  below  the  HfSiON,  or  at  the  interface  between  the  SiC>2  and  high-K  layer.  We  further 
stipulate  that  the  positive  charge  in  that  region  is  highly  mobile  in  the  lateral  sense.  This 
hypothesis  is  not  inconsistent  with  lateral  charge  redistribution  effects  observed  in  other 
dielectric  system  such  as  in  the  SONOS  memory  cell  [10,  11].  The  mobile  charge  may  be  due  to 
holes,  or  perhaps  protons;  we  will  assume  here  that  the  charge  is  due  to  holes.  We  suggest  that 
the  drain  voltage  produces  a  lateral  field  in  both  the  semiconductor  and  in  the  dielectric  layers 
that  removes  the  holes  by  conducting  them  out  via  the  drain,  thus  reducing  the  hole  population  in 
the  dielectric  stack  by  an  amount  that  increases  with  increasing  VdS.  Because  the  drain  voltage  is 
not  “felt”  far  from  the  drain,  the  mechanism  must  involve  removal  of  holes  near  the  drain, 
followed  by  transport  of  holes  toward  the  drain  from  farther  inside  the  channel,  much  as  charge 
is  conducted  out  of  the  drain  from  the  Si  inversion  layer  during  normal  device  operation. 


3.4.  Estimation  of  the  magnitude  of  different  contributions  to  AVy, 
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The  simplest  approximation  to  model  the  NBTI  induced  charging  is  to  assume  two 
distinct  “species”,  interface  states  and  charges  trapped  in  the  gate  dielectric  [8].  With  this 
assumption  we  propose  a  simple  qualitative  model  of  NBTI  charging  kinetics  for  the  devices  we 
have  studied  that  actually  allows  us  to  extract  the  time  dependent  variation  of  the  two 
components  of  AVti,(NBTI).  We  begin  by  noting  that  if  we  take  the  data  of  Fig.  1  and  plot  AVth 
at  fixed  time  as  a  function  of  VdS,  we  can  generate  a  series  of  curves  as  shown  in  Fig.  4a,  which 
emphasize  the  variation  of  AVt|,(NBTI)  with  Vds.  One  observes  that  the  curves  are  in  fact  not 
parallel,  but  tend  to  diverge  at  lower  VdS.  We  now  assume  that  (i)  holes  are  trapped  rapidly  in 
comparison  with  the  generation  of  interface  states,  so  that  for  early  times  the  interface  state 
component  is  negligible,  and  (ii)  that  interface  state  generation  is  roughly  independent  of  drain 
voltage.  One  can  now  write  the  trapped  hole  and  interface  state  contribution  to  the  total  threshold 
voltage  shift  as 

AVtll(0  =  AFchh(O  +  AF“(0.  (6) 

It  is  reasonable  also  to  suggest  that  for  times  tj,  t2  (t2  >  ti)  not  too  different  from  one  another, 
hole  trapping  increases  proportionally  with  time  so  that 

AFfhft2)=  a(OAltt<tt)  +  .  (7) 

This  equation  models  the  divergence  of  the  curves  at  lower  VdS:  for  longer  times,  more  interface 

states  are  generated.  In  Eq.  (7),  (t3)  Al'/h  (t3)  is  the  shift  due  to  interface  states  generated 

between  times  t2  and  t|.  Finally,  if  we  take  the  curve  in  Fig.  4a  corresponding  to  a  time  of  20  s  as 
one  for  which  ^'th (an  arbitrary  choice,  provided  t  is  short),  we  can  construct  the  interface 
state  density  as  a  function  of  time.  Plots  of  AFth(ts)AViii(t2  )  as  a  function  of  rVft)Avn>(ti)  for 
two  choices  of  t2  and  t|  yield  straight  lines,  as  shown  in  Fig.  4b,  with  slope  a(t)  and  intercept 
A*  rrfcfts)A1/tss (f).  The  density  of  interface  states  determined  from  the  intercepts  is  plotted  as  a 
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function  of  time  in  Fig.  5.  As  mentioned  previously,  implicit  in  this  calculation  is  that  interface 
state  generation  is  independent  of  VdS,  which  is  reasonable  given  that  the  field  (and  G(E)  from  B. 
above)  is  uniform  across  the  large  fraction  of  the  channel.  In  choosing  ft  and  t2  values,  we  found 
that  a  difference  of  no  more  than  1 00  s  gives  reasonably  straight  line  plots  such  as  those  in  Fig. 
4b;  accordingly  we  chose  t  =  20s,  50s,  100  s,  200  S...600  s.  Given  the  level  of  approximation 
adopted  here,  we  feel  that  the  results  are  reasonable  and  that  the  model,  although  qualitative,  is 
justified. 

In  developing  the  plot  of  interface  states  shown  in  Fig.  5  we  have  assumed  certain  things 
about  NBT1  charging:  that  it  comprises  separate  hole  and  interface  state  generation  in  HfSiON, 
and  that  the  interface  states  are  generated  relatively  slowly,  and  after  holes  are  trapped.  The  first 
of  these  implications  is  consistent  with  the  literature,  although  in  fact  it  is  more  likely  that  a  wide 
range  of  interface  and  near-interface  defects,  with  widely  varying  response  times,  exists.  In  that 
case,  what  we  (and  others  who  write  equations  similar  to  Eq.  (6))  have  done  is  to  probe  two 

broad  classes  of  defect:  one  that  is  relatively  slow  (AE A )  and  one  that  is  relatively  fast  (hVrl£  y 
The  latter  assumption  we  cannot  prove,  but  we  note  that  it  is  consistent  with  interface  state 
generation  models  in  which  H2  diffusing  through  the  oxide  “cracks”  on  a  positively  charged 
center  and  de-passivates  a  silicon  dangling  bond  [12].  That  model  has  been  developed  to  explain 
interface  state  generation  in  radiation  damaged  oxides,  and  we  suggest  that  it  may  be  appropriate 
here  as  well. 


4.  Conclusion 

We  have  performed  NBTI  measurements  in  high-k  gate  stack  devices,  and  these  reveal  an 
anomalous  dependence  of  the  threshold  voltage  shift  on  Vjs  as  one  changes  the  stress  conditions 
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from  the  linear  to  the  saturation  mode.This  effect  was  not  observed  in  devices  with  pure  SiC>2 
gate  dielectrics.  Customarily,  NBTI  in  MOSFETs  is  studied  either  with  the  source,  body  and 
drain  contacts  shorted  or  with  a  small  potential  (~  0.1  V  or  less)  between  the  source  and  drain 
contacts.  We  are  presently  unable  to  say,  on  the  basis  of  our  measurements  and  of  published 
data,  whether  or  not  the  observed  effect  is  particular  to  our  devices  or  generic  to  other  devices  or 
dielectric  gate  stacks.  If,  for  example,  the  effect  is  generic  to  high-k  gate  stacks  involving 
complex  oxides,  our  data  has  important  implications  for  estimates  of  the  reliability  lifetime  due 
to  NBTI  in  p  channel  MOSFETs.  We  provide  evidence  that  this  effect  is  consistent  with  removal 
of  holes/protons  at  the  drain  contact  by  conduction  along  the  dielectric  stack,  i.e.  lateral  charge 
transport. 

A  simple  two-species  model  that  assumes  interface  state  buildup  is  slow  compared  to 
hole  trapping  allows  us  to  extract  values  for  the  interface  state  density  component  of  AV(|,(NBTI) 
in  the  HfSiON  based  devices.  The  interface  state  density  tends  to  saturation  at  times  greater  than 

it 

500  s,  at  which  time  it  is  found  to  be  on  the  order  of  1 0  cm'  .  These  values  are  consistent  with 
published  data  on  interface  state  generation  by  NBTI  in  devices  having  a  variety  of  dielectric 
gate  stacks  [9,  13,  14].  Furthermore,  they  are  consistent  with  the  values  one  can  deduce  from  the 
data  in  Fig.  2  for  pure  SiC>2  gate  devices  where  we  expect  interface  state  generation  by  NBTI  to 
dominate  charge  trapping  in  the  “bulk”  of  the  dielectric. 

Further  studies  are  clearly  relevant  to  ascertain  how  widespread  the  phenomena  is  that  we 


have  observed  in  our  devices 
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Figure  Captions 

Fig  L  a)  Normalized  change  in  drain  current  in  HfSiON  based  devices  measured  for  NBTI  stress 
at  Vgs  =  -2.0  V,  T  =  175  °C,  and  varying  VdS  as  indicated.  Data  was  taken  approximately  every 
0.75  s  but  for  clarity  only  every  10lh  point  is  displayed.  Measurements  were  made  at  VdS  =  -1.5 
V  as  well  but  for  clarity  are  not  shown,  b)  Threshold  voltage  shifts  determined  from  the  data  of 
Figure  la  using  Eq.  (3). 

Fig.  2.  a)  Alds(t)  /Ids°  for  SiC>2  gate  p-channel  MOSFETs  at  195  °C.  Measurements  were  made  in 
the  linear  regime  (  □)  and  saturation  regime  (O)  b)  Data  converted  to  AVn,(t)  using  Eq.  3.  In 
both  cases  the  gate-source  voltage  was  maintained  at  -2V. 

Fig.  3.  a)  Vertical  electric  field  Ez(y)  in  the  oxide  as  a  function  of  distance  along  the  channel 
obtained  from  a  quasi-2D  model  of  the  surface  potential.  Solid  line:  VdS  =  -  0.1  V;  dashed  line: 

Vds  =  -  2.0  V.  b)  Generation  rate  fj(Ez)~Eze)  G(E)~EeyE  in  arbitrary  units  for  VdS  =  -  2.0  V. 

Fig.  4.  a)  AVth  in  HfSiON  devices  plotted  at  selected  times,  as  indicated,  and  as  a  function  of  the 
drain  voltage  VdS  applied  during  stress.  The  lines  are  guides  for  the  eye.  (b)  AV,h(t2)  as  a  function 
of  AV,i,(t|),  for  three  sets  of  t2,  tj.  The  intercept  along  the  AVth(t2)  axis  in  each  case  is  due  to 
interface  states  generated  between  times  fi  and  t2.  At  time  fi  =  20s,  interface  state  density  is 
assumed  to  be  0. 

Fig.5.  Density  of  interface  states  determined  from  the  intercepts  of  curves  such  as  those  in  Figure 
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